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We investigate resonant capture of small bodies by planets that migrate inwards, 
D . using analytic arguments and three-body integrations. If the orbits of the planet 

in 



and the small body are initially circular and coplanar, the small body is captured 
when it crosses the 2:1 resonance with the planet. As the orbit shrinks it becomes 
more eccentric, until by the time its semimajor axis has shrunk by a factor of 
four, its eccentricity reaches nearly unity (1 — e ^C 10~ 4 ). In typical planetary 
systems, bodies in this high-eccentricity phase are likely to be consumed by the 
CSj , central star. If they can avoid this fate, as migration continues the inclination 

o. 

axis is a factor of nine smaller than at capture, at which point the small body 
is released from the 2:1 resonance on a nearly circular retrograde orbit. Small 
bodies captured into resonance from initially inclined or eccentric orbits can also 
be ejected from the system, or released from the resonance on highly eccentric 
polar orbits (i ~ 90°) that are stabilized by a secular resonance. We conclude 
that migration can drive much of the inner planetesimal disk into the star, and 
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that post- migration multi-planet systems may not be coplanar. 

Subject headings: celestial mechanics, stellar dynamics — solar system: formation 

1. Introduction 

Giant planets are found orbiting nearby stars with semimajor axes as small as 0.04 AU 
( [Schneider 2000] ). The difficulty of forming planets so close to a star has led to the suggestion 
that these planets formed at larger radii and migrated inwards, perhaps through gravitational 
interactions with the gaseous protoplanetary disk or planetesimal disk (|Lin, Bodenheimer fc 
Richardson T996J [Murray et al. 1998Q . 



When massive planets migrate, they can capture smaller bodies in mean-motion reso- 
nances. Resonant capture in the solar system was invoked by Goldreich (1965) to explain 
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the frequent occurrence of commensurable periods among planetary satellites; in this case 
the satellites are migrating outwards due to tidal friction (see [Peale 1999| for a review). 



Similarly, the outward migration of Neptune probably led to the capture of Pluto and many 
Kuiper-belt objects ("Plutinos") into the 2:3 resonance with Neptune (Malhotra 1993, 1995). 

All known examples of resonant capture in the solar system involve outward migration. 
However, the compact orbits of exosolar planets are believed to result from inward migration. 
The purpose of this paper is to explore some of the novel features of resonant capture during 
inward migration. 

We shall use the Delaunay elements, 

L = (M„a) 1 / 2 , G = L(l-e 2 ) 1/2 , H = Gcosi, £, g = u, h = Q; (1) 

here M* is the mass of the star in units where the gravitational constant is unity, a, e, and 
% are the semimajor axis, eccentricity, and inclination; £, uo and Q are the mean anomaly, 
argument of periastron, and longitude of the ascending node; G is the specific angular 
momentum and H is the z-component of the angular momentum. In these variables the 

Kepler Hamiltonian is 

M 2 
Hk(L) = -^, (2) 

and the mean motion n = 8H K /dL = M^/L 3 = (M^/a 3 ) 1 / 2 . 

For simplicity we shall consider only the case of a test particle that is captured by a 
massive body on a circular orbit with zero inclination. We shall call this body "Jupiter" ; its 
mass, semimajor axis, and mean motion are denoted by Mj, aj, and nj. 



2. Analytic results 

A (p, q) mean-motion resonance occurs when pn ~ (p + q)nj, where p > and q > —p are 
integers. We distinguish inner resonances, which lie inside Jupiter (a < aj, n > nj, q > 0), 
from outer resonances (a > aj, n < nj, q < 0). 

The Hamiltonian of the test particle can be written in the form 

H K (L) + Y, A pgm (L, G, H) cos[p£ + (p + q)(g + h - <f>j) - 2mg}; (3) 

p,q,m 

here <f>j(t) = fnj(t)dt is the azimuth of Jupiter, and p, q, and m are integers. The angles h 
and (j)j can only appear in the combination h — <pj because the Hamiltonian is invariant to 
changes in the origin of the azimuth. The restriction that the coefficient of g must be an even 
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integer 2m arises because the planetary potential is symmetric about the equatorial plane, so 
that the Hamiltonian must be invariant under the transformation /i— >/i + 7r, g— >g + 7r. For 
low-eccentricity, low-inclination orbits, |A pgm | oc (Mj/M*)e'' 3_2m 'i' 2m ' forp, q ^ 0. The terms 
with p + q = can also contain a contribution from an additional axisymmetric potential, 
due for example to a protoplanetary disk or the equatorial bulge of the central star. Terms 
with p = q = are sometimes called the secular Hamiltonian since they do not depend on 
either the mean anomaly £ or the azimuth of Jupiter <pj and thus are time-independent for 
Keplerian motion. 

Each (p, q) mean-motion resonance is represented in the Hamiltonian (^) by an infinite 
series in m. We can isolate the slowly varying resonant angle by a canonical transformation 
to new momenta and coordinates (R, r) defined by the mixed-variable generating function 

S(R, £, g, h) = Rd + R 2 g + R 3 [p£ + (p + q )( g + h- 0j)]. (4) 

We have 

dS dS dS 

L = — = R 1 +pR 3l G = — = R 2 + (p + q)R 3l H = — = (p + q)R 3l (5) 

OS dS dS 

n = 7^- = ^, r2 = ~Q R ~ = g i r3 = 'q r ~ =K + (P + <?)(# + ^-0j)- (6) 

Near the resonance, r\ varies much more rapidly than r 2 or r 3 , so we can average the 
Hamiltonian over the fast variable r\ at fixed r 2 , r 3 . The conjugate action R\ is a constant 
of motion in the averaged Hamiltonian and an adiabatic invariant of the original Hamiltonian. 
Therefore 

C„ = (p + q)Ri = (p + q)L-pH = (M*a) l l 2 [{p + q) - p(l - e 2 ) 1 ' 2 cosi] (7) 

is adiabatically invariant. 

The averaged Hamiltonian has two degrees of freedom. The nature of the trajectories 
in the averaged Hamiltonian depends on the relative strength of the secular and resonant 
Hamiltonians. If the secular Hamiltonian is either large or small compared to the resonant 
Hamiltonian, the timescales for libration or circulation of the associated coordinates r 2 and 
r 3 are well-separated, and the problem can be reduced to one degree of freedom by a second 
averaging. In the intermediate case, resonances with different m can overlap, leading to 
chaotic motion. All three regimes are present in various solar-system contexts. 

Despite these complexities, some of the principal features of resonant capture can be 
deduced directly from the adiabatic invariant C pq (eq. 0). Suppose that the test particle is 
initially on a circular, equatorial orbit at semimajor axis aj. Then C pq = (M+diY^q. Once 
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the particle is captured into resonance, its semimajor axis will migrate along with Jupiter's 
to preserve the mean-motion resonance, thereby pumping up the test particle's eccentricity 
and/or inclination. The adiabatic invariance of C pq then requires that 

k = (1 - e 2 ) 1/2 cos* = 1 + 2 [l - ( ai /a) 1/2 ] . (8) 

The quantity on the left is initially unity (since e = i = 0) and can only decrease. For an 
outer resonance (q < 0) this requires that a t /a < 1; in other words capture into an outer 
resonance can only occur if the perturber migrates outward. If the perturber migrates a long 
way outwards, then a^ja — > and k — > /too = 1 + q/p- Note that < Kqo < 1 since p > 0, 
— p < q < 0; thus captured orbits never achieve radial (e = 1) or polar (i = ^tc) orbits. 

Similarly, capture into an inner resonance (q > 0) can only occur if the perturber 
migrates inward. However, in this case the eccentricity and inclination evolution is more 
dramatic. The quantity k — > when a = q 2 cii/(p + q) 2 ; in other words, by the time Jupiter 
has migrated inward by a factor of q 2 /(p + q) 2 after capture, the test particle has been 
driven onto either a radial or a polar orbit. If migration continues past this point and the 
test particle remains in resonance, it must be driven onto a retrograde orbit (« < 0). The 
minimum possible value of k is —1 (for e = 0, i = tt). This is attained when Jupiter has 
migrated inward by a factor of q 2 /(2p + q) 2 since capture; if migration continues past this 
point the test particle cannot remain in resonance. 

These simple considerations suggest that inward migration of a planet can excite small 
planets or planetesimals into high-eccentricity, high-inclination, or retrograde orbits. In 
the remainder of this paper we shall investigate the evolution and fate of these objects 
numerically. 



3. Numerical integrations 

In our integrations we use units in which the gravitational constant, Jupiter's initial semima- 
jor axis a,ji and orbital period are all unity, so M* = 4ti 2 . Jupiter's mass is Mj = 0.001M*. 
Jupiter is assumed to follow a circular orbit that migrates inward according to the rule 

a,j(t) = ajf + (1 - a Jf ) exp(-£/V), (9) 

where ajf = 0.1 and r = 1000. The test-particle integrations were carried out in regularized 
coordinates. 

Figure [I] shows the evolution of a test particle initially on a zero-inclination, circular 
orbit with semimajor axis a = 0.59. The test particle is captured into the 2:1 resonance 
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Fig. 1.— Resonant capture of a test particle by an inward migrating planet. From top 
to bottom, the panels show the ratio of mean motions of Jupiter and the test particle, the 
semimajor axis (two plots on different scales), the eccentricity, and the inclination of the test 
particle. The test particle is initially on a circular, equatorial orbit (e« = 0, %i — 0), shown 
as an enlarged open circle in Figure |2|. It is captured into the 2:1 resonance with Jupiter and 
later released on a retrograde near-circular orbit (i = 180°, e = 0.09). 
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(p = q = 1) with Jupiter at t/r = 0.073. Thereafter its orbit shrinks and its eccentricity 
grows, as predicted by the adiabatic invariant k (eq. ||). Eventually at t/r = 1.90, when the 
semimajor axis has shrunk by a factor of q 2 /(p + q) 2 = 4 since capture, the orbit is nearly 
radial (1 — e <C 10~ 4 ). In typical planetary systems, a small body in this high-eccentricity 
phase would be vaporized by or collide with the star, but we shall continue to follow the 
dynamical evolution assuming that the particle survives. At this point the inclination flips 
rapidly from to 7r (recall that inclination is undefined for radial orbits). Once the inclination 
is 7r the eccentricity shrinks as the orbit continues to shrink, again as predicted by adiabatic 
invariance. Finally, at t/r = 5.40, when the semimajor axis is a factor of (1 + 2p/q) 2 = 9 
smaller than at capture, the test particle is released from the resonance on a retrograde, 
equatorial, near-circular orbit. 

Figure [2] shows the final states of an ensemble of test particles, on initially circular, zero- 
inclination orbits with initial semimajor axes Oj between 0.26 and 0.6. The planet migrates 
from aji = 1 to ajf = 0.1, so the 2:1 resonance migrates from am = 0.6300 to or/ = 0.0630. 
According to equation (|8|), particles will be captured into resonance during the migration if 
a i > a>Rf, flipped onto retrograde orbits if a, t > (1 + p/q) 2 aRf = 4or/ = 0.2520, and released 
from resonance on nearly circular, retrograde orbits if Oj > (2p/q+ l) 2 a#/ = 9aRf = 0.5670. 
For particles with 4auf < a« < 9aRf the eccentricity is given by 



iRf 



:i-4) 1/2 



2 



(10) 



which is plotted as a dotted line. All of these predictions agree very well with the results 
shown in the Figure. 

The evolution is more complicated when the test particles have non-zero eccentricities 
and inclinations. Figure |3] shows the final states of test particles with the same range of 
initial semimajor axes as in Figure 0, but with eccentricities and inclinations chosen from 
the Schwarzschild or Rayleigh distribution, 

( e 2 i 2 \ 
n(e,i)dedi oc exp r ^ edeidi. (11) 

V 2e o 2 *o/ 
At the end of the integration, most of the particles are found in one of three final states: 

• The particles shown as filled circles are in retrograde (i ~ n) orbits trapped in the 2:1 
resonance, like those shown in Figure |2|. If migration had continued, these would even- 
tually be released from resonance on nearly circular orbits, after their semimajor axis 
had migrated by a factor 9 since capture (this has already happened for two particles, 
marked by open circles). These particles typically have small initial eccentricities and 
inclinations, e^ <; 0.05, i, ^ 5°. 



-7- 



150 

w 
0) 
0) 
u 100 

■~ 50 




1 

0.8 

0.6 

0.4 

0.2 





Q) 



• •• • • • 



-i i i |- 
■ •«••• 



• © 



0^~ 



e^O, i,=0 



-•»•••- 



j_ 



j_ 



© 



0.3 



0.4 0.5 

a. 



0.6 



Fig. 2. — Resonant capture of test particles on initially circular, equatorial orbits. The top 
and bottom panels show the final inclination and eccentricity, if and e/, as a function of 
the initial semimajor axis o^. Filled circles represent test particles that are still in the 2:1 
resonance at the end of the integration, while open circles represent particles that have been 
released from resonance. The dotted line shows the prediction of equation ([H]) for the final 
eccentricity. The enlarged open circle is the test particle whose evolution is plotted in Figure 
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(b) e =0.06, i =0 
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Fig. 3.— Evolution of test particles from low-eccentricity, low-inclination orbits. The 
initial eccentricities and inclinations are chosen from the Schwarzschild distribution (|TT|), 
with characteristic widths e = 0.03, i = (top panels), e = 0.06, i = (middle panels), 
and eo = 0, i$ = 5° (bottom panels). The left- and right-hand columns show the final 
inclination and eccentricity as a function of initial eccentricity or inclination and initial 
semimajor axis. The dotted lines in the right-hand panels show the prediction of equation 
(|10|) . The symbols denote the final states of the test particles at the end of the integration: 



circles represent retrograde orbits 
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triangles represent nearly polar orbits 



vi- 



and crosses represent escape orbits or chaotic orbits that are likely to escape in the future. 
Filled symbols denote particles that are in the 2:1 resonance at the end of the integration. 
Eccentricities and inclinations are computed in the star-centered frame, except for escape 
orbits, where the elements are computed in the center-of-mass frame. The enlarged symbols 
in the bottom panels are the test particles whose evolution is plotted in Figures [| and ||. 
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The particles marked by crosses are on escape orbits, or chaotic orbits that are likely to 
escape in the future. In most cases these particles are pumped onto high-eccentricity 
orbits, released from resonance, and then suffer a close encounter with Jupiter, leading 
to a random walk in semimajor axis and eventual escape. The evolution of the orbital 
elements in a typical case is shown in Figure |j. These particles typically have initial 
eccentricities e« ^ 0.05. Note that the final eccentricity and inclination for the escaping 
orbits are given in the center-of-mass frame, while all other orbital elements in this and 
other Figures are given in the star-centered frame. 

The particles marked by open triangles were trapped in the 2:1 resonance and later 
released on nearly polar orbits (i ~ |7r). These particles have high eccentricities and 
arguments of perihelion u that librate around ±^vr. This secular resonance protects 
high-eccentricity orbits from close encounters with Jupiter, so the particle orbits are 
stable on the timescale of our integrations (1 x 10 4 initial orbital periods of Jupiter). 
These orbits typically have initial inclination z'j ^ 5°. An example is shown in Figure 



4. Discussion 

Giant planets are found orbiting close to many nearby stars. Because the conditions in 
the inner parts of protoplanetary disks are unfavorable to planet formation, it is likely that 
these planets formed at much larger radii and migrated inward. Moreover, because the 
characteristic evolution times in a protoplanetary disk are longer at larger radius, it is likely 
that planetesimals and small planets had already formed interior to this giant planet before 
it began its inward journey. In these circumstances resonant capture by the inward migrating 
planet is likely to occur at the 2:1 mean-motion resonance. 

We have shown that small bodies that are captured from circular, equatorial orbits into 
the 2:1 resonance will be excited onto nearly radial orbits by the time their semimajor axis 
has shrunk by a factor of four, and released from the resonance on nearly circular retrograde 
orbits when the semimajor axis has shrunk by a factor of nine. In most of the known 
planetary systems, these small bodies would be consumed by the central star during this 
high-eccentricity phase. Thus, substantial migration of a giant planet is likely to pollute the 
outer parts of the star with most of the mass in the interior planetesimal disk: specifically, 
if the planet migrates from aji to ajf < \cLji, then the planetesimals with initial semimajor 
axes in the range 2. 520a j/ < a < 0.630ajj will be excited to high-eccentricity orbits and are 
likely to be vaporized by or collide with the central star. 
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Fig. 4.— The evolution of a test particle with initial eccentricity and inclination e« = 0, 
ii = 5.45°, marked by an enlarged cross in Figure || From top to bottom, the panels show 
the ratio of mean motions of Jupiter and the test particle, the semimajor axis (two plots 
on different scales), the eccentricity, and the inclination of the test particle. The particle 
is temporarily captured into the 2:1 resonance and later released on a chaotic orbit, which 
eventually escapes after a close encounter with Jupiter. The eccentricity and inclination are 
not constant at large time because we are working in a frame centered on the star rather 
than the center of mass of the star and Jupiter. The final eccentricity and inclination in the 
center-of-mass frame are e = 1.01, i = 69°. 
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Fig. 5.— The evolution of a test particle with initial eccentricity and inclination e = 0, 
i = 2.37°, marked by an enlarged triangle in Figure |3|. From top to bottom, the panels show 
the ratio of mean motions of Jupiter and the test particle, the semimajor axis (two plots 
on different scales), the eccentricity, the inclination, and the argument of periastron (the 
angle in the orbital plane from the equatorial plane to the periastron). The test particle is 
captured into the 2:1 resonance and later released on a polar (i ~ ~7r) orbit. The particle 
is protected from close encounters with the planet by a secular resonance, in which the 
argument of periastron librates around \ r n. 
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Bodies that do not collide with the star will eventually be released from resonance as 
migration continues. Some of these will be ejected, but some are protected from ejection by 
secular resonance and can survive on polar, high-eccentricity orbits. Others are released from 
resonance on circular and equatorial but retrograde orbits, and will survive if the migration 
stops before the giant planet reaches their semimajor axis. Thus, if the captured planets 
survive the high-eccentricity phase, their post-migration configuration may not be planar; 
if the orbits are coplanar, they may not all orbit in the same direction. If multiple bodies 
formed from a single disk ("planets") are not necessarily on coplanar orbits, they may be 
difficult to distinguish from bodies formed from fragmentation of a collapsing gas cloud 
("stars"). Radial- velocity observations cannot directly measure the relative inclinations of 
planets in multi-planet systems; this determination requires astrometric measurements such 
as those planned for the Space Interferometry Mission (SIM). 

This preliminary study has not addressed many important issues, including the effects 
of an eccentric planet orbit, the stability of these resonances to gravitational noise from 
other planets and massive planetesimals, or migration in systems containing two or more 
giant planets. 

This research was supported in part by NASA Grant NAG5-7310. 
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